#include <stdio.h>
#include <math.h>
float a[50][50], b[50][50],c[50][50];

float determinante(float[50][50]);

float determinante(float matriz[50][50])
{
	int i,j,k,l;
	det=matriz[1][1];
	for(k=1; k<n; k++) 
	{
		l=k+1;
		for(i=l;i<=n;i++)
			for(j=l;j<=n;j++)
				matriz[i][j] = ( matriz[k][k]*matriz[i][j]-matriz[k][j]*matriz[i][k] )/matriz[k][k];
		det = det*matriz[k+1][k+1];
	}
	return det;
}